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Leptogenesis through direct infiaton decay to light particles 
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We present a scenario of nonthermal leptogenesis following supersymmetric hybrid inflation, in the 
case where infiaton decay to both heavy right handed neutrino and SU(2)l triplet superfields is 
kinematically disallowed. Lepton asymmetry is generated through the decay of the infiaton into 
light particles by the interference of one-loop diagrams with right handed neutrino and SU(2)l 
triplet exchange respectively. We require superpotential couplings explicitly violating a U(l) R- 
symmetry and R-parity. However, the broken R-parity need not have currently observable low- 
energy signatures. Also, the lightest sparticle can be stable. Some R-parity violating slepton decays 
may, though, be detectable in the future colliders. We take into account the constraints from 
neutrino masses and mixing and the preservation of the primordial lepton asymmetry. 

PACS numbers: 98.80.Cq, 12.10.Dm, 12.60.Jv 



I. INTRODUCTION 

The standard model (SM) of electroweak and strong in- 
teractions, despite its many successes, leaves unanswered 
many questions in particle physics and cosmology. In 
particular, it does not address the following problems: 
the origin of electroweak symmetry breaking and of the 
hierarchy between the electroweak scale M\y and the (re- 
duced) Planck scale mp ~ 2.44 x 10 18 GeV; the origin 
and size of neutrino (y) masses and mixing; the cosmolog- 
ical horizon and flatness problems and the origin of den- 
sity perturbations; the generation of the observed baryon 
asymmetry of the universe (BAU) . 

The problem of stabilizing the electroweak scale rela- 
tive to the fundamental scale Ms (which we take to be 
close to the Planck scale) is solved by softly-broken su- 
persymmetry (SUSY). For definiteness, we take the case 
of gravity-mediated SUSY-breaking, where the gravitino 
mass is of order 1 TeV. In this framework, the vac- 
uum expectation values (VEVs) of the electroweak Higgs 
superfields h%, h% can be determined by radiative cor- 
rections to the scalar potential 0. However, this ex- 
planation of the electroweak symmetry-breaking requires 
a term \ih\hi in the superpotential with [i ~ 1 TeV, 
whereas a priori the value of /i is expected to be of 
the order of the fundamental scale. This fi problem of 
the minimal SUSY standard model (MSSM) can be ad- 
dressed by imposing a symmetry that forbids the above 
fj, term, which is then broken in a controlled fashion. A 
R-symmetry, broken by the soft SUSY-breaking terms, 
can fulfill this role (see below). 

The smallness of neutrino masses is naturally explained 
in a model with additional SM singlet chiral fermions 
(v c ), the right handed neutrinos (RHNs), since these can 
obtain heavy Majorana masses (~ M u e) and also par- 
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ticipate in Yukawa couplings with the SU(2)l doublet 
neutrinos resulting in Dirac neutrino masses (~ Tod) 
after the electroweak symmetry breaking. The result- 
ing seesaw mass matrix [jjj has small eigenvalues of or- 
der mr D /M v c, which can be of the right order of magni- 
tude to explain atmospheric and solar neutrino oscilla- 
tions, while being consistent with cosmological bounds 
on neutrino masses. Another mechanism of inducing 
nonzero neutrino masses is the introduction of heavy 
scalar SU(2)l triplets T with lepton number —2, coupling 
to the SU(2)l doublet lepton fields. After the breaking 
of the electroweak symmetry, the neutral components of 
these triplets can acquire small VEVs inducing Majorana 
neutrino masses . In general, these two mechanisms for 
generating neutrino masses can coexist. 

The cosmological horizon and flatness problems are 
most elegantly solved by inflation, which can also gen- 
erate the primordial density perturbations required for 
structure formation in the universe p|. Moreover, infla- 
tion, which can be easily incorporated in realistic particle 
physics models, is strongly favored by the recent data 
on the angular power spectrum of the cosmic microwave 
background radiation (CMBR). 

If neutrinos get mass either by coupling to heavy SM 
singlet fermions or to heavy scalar SU(2)/, triplets, it may 
be possible to generate 0, a primordial lepton asym- 
metry in the out-of-equilibrium decay of these heavy de- 
grees of freedom, in case they were thermally produced in 
the early universe. This asymmetry is then reprocessed 
at the electroweak phase transition to yield the observed 
BAU. However, there is a tension between correct neu- 
trino masses and this thermal leptogenesis scenario in 
SUSY models because of the gravitino problem 0, 0] . 
With a gravitino mass of around 1 TeV, the reheat tem- 
perature T r h should not exceed 10 9 GeV since graviti- 
nos produced in thermal scattering processes would de- 
cay late presumably into photons and photinos and, if 
sufficiently numerous, interfere with the successful pre- 
dictions of standard big bang nucleosynthesis. On the 
other hand, one also requires that the heavy fields whose 
decay creates lepton asymmetry be generated in sufficient 
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numbers. So, their masses should not exceed T r h, which 
leads to unacceptably large light neutrino masses. How- 
ever, this problem can be alleviated by assuming 
that there is some degree of degeneracy between the rele- 
vant RHNs or heavy SU(2)l triplets, which enhances the 
generated lepton asymmetry, and perhaps also that the 
branching ratio of the gravitino decay into photons and 
photinos is less than unity, which somewhat relaxes fl(i| 
the gravitino constraint on T r h- 

The tension between v masses and the gravitino prob- 
lem can be more naturally relaxed by considering non- 
thermal leptogenesis |l3[ at reheating. However, in exist- 
ing scenarios where the inflaton decays into RHN 
or SU(2)i triplet superfields, this still puts a restriction 
on the masses of these particles: the decay products of 
the inflaton must be lighter than half its mass m m { . Lep- 
ton asymmetry is generated in the subsequent decay of 
the RHN or SU(2) L triplet superfields. 

In this work, however, we consider the consequences 
of allowing all the RHN and/or SU(2)l triplet super- 
fields of the model to be heavier than mi n f/2 (see also 
Ref. 0). Lepto genesis could then occur only via the 
direct decay of the inflaton to light particles (sec also 
Ref. 0). We take a simple SUSY grand unified the- 
ory (GUT) model naturally incorporating the standard 
SUSY realization [3 [13 of hybrid inflation [U, which 
does not require tiny parameters and is, undoubtedly, one 
of the most promising inflationary scenarios. (For exten- 
sions of standard SUSY hybrid inflation, see Ref. |2l|-) 
In global SUSY, the flatness of the inflationary path at 
tree level is guaranteed by a U(l) R-symmetry. The r\ 
problem of sizable supergravity (SUGRA) contri- 
butions to the inflaton mass on the inflationary path, 
which could easily invalidate inflation, is reduced, in this 
case, to controlling the magnitude of a single term in the 
Kahler potential [22] (see also Ref. [23j). Finally, radia- 
tive corrections provide |19| a logarithmic slope along the 
classically flat direction, needed for driving the inflaton 
towards the SUSY vacua. 

The uproblem is solved by employing the mechanism 
of Ref. 0|. (For an alternative solution of the /i prob- 
lem, see Ref. 01 •) The global R-symmetry of the model 
forbids the appearance of a /i term in the superpoten- 
tial. On the contrary, it allows the existence of the tri- 
linear term Sh\h2, where S is the gauge singlet infla- 
ton of the standard SUSY hybrid inflation. After the 
GUT gauge symmetry breaking, the soft SUSY-breaking 
terms, which generally violate the R-symmetry, give rise 
to a suppressed linear term in S and, thus, this field ac- 
quires a VEV of order the electroweak scale divided by 
a small coupling constant. The above trilinear coupling 
can then yield a /i term of the right magnitude. 

The inflaton consists of two complex scalar fields with 
tree-level couplings to the electroweak Higgses and Hig- 
gsinos derived from the above trilinear term. After the 
end of inflation, it oscillates about the SUSY vacuum and 
eventually decays predominantly into electroweak Higgs 
superfields via these tree-level couplings, thereby reheat- 



ing the universe. We find that, in the case of one fermion 
family, both heavy RHN and SU(2)l triplet superfields 
are necessary if diagrams producing a nonzero lepton 
asymmetry are to exist. Since these heavy fields can 
only appear in intermediate states of the inflaton decay, 
we must create the asymmetry directly from this decay. 
Indeed, leptogenesis can occur in the subdominant decay 
of the inflaton into lepton and Higgs superfields through 
the interference between different one-loop diagrams with 
RHN and SU(2) l triplet exchange respectively. The lep- 
ton asymmetry is proportional to a novel CP-violating 
invariant product of coupling constants. 

For a nonzero asymmetry, we also need to include 
some couplings in the superpotential that explicitly vio- 
late both the U(l) R-symmetry and its Zi matter parity 
subgroup, which remains unbroken by the soft SUSY- 
breaking terms. Although these couplings involve super- 
heavy fields, the matter parity violation may have some 
observable consequences at low energy such as the pos- 
sible instability of the lightest SUSY particle (LSP). In- 
deed, if this particle contains a Higgsino component, it 
could decay predominantly into a pair of Higgses and 
a lepton. This channel can, though, be easily blocked 
kinematically if the LSP is not too heavy (see Sec. IVj l. 
On the contrary, some R-parity violating slepton decays 
which may be detectable in the future colliders are typi- 
cally present (see Sec. IVIII) . 

We find that the value of the BAU from the Wilkinson 
microwave anisotropy probe (WMAP) data |6| can be 
easily achieved given constraints from other observables, 
notably the reheat temperature and neutrino masses and 
mixing, and CP- violating phases of order unity. However, 
the constraint from v masses and mixing is considerably 
weakened by the existence of two separate contributions, 
namely those of the usual seesaw mechanism and from 
the Higgs SU(2) L triplets (see alsoc Ref. On the 

contrary, the requirement that the initial lepton asym- 
metry is protected from lepton number violating 2^2 
scattering processes which are in equilibrium in the early 
universe imposes very stringent constraints on the pa- 
rameters of the model, which we also take into account 
in our analysis. The prediction for the spectral index of 
density perturbations is typical of SUSY hybrid inflation 
models (see e.g. Ref. j2j|). 

Thus, an acceptable value of the BAU can be ob- 
tained within a consistent model of cosmology and parti- 
cle physics, without requiring additional fine-tuned cou- 
pling constants and without necessarily putting strong 
constraints on observables such as neutrino masses and 
mixing. Moreover, although the scenario requires viola- 
tion of the R-symmetry, it is not necessary to introduce 
superpotential terms which would lead to currently ob- 
servable R-symmetry violating effects. 

In Sec. [n] we introduce our model and describe some 
of its salient features. In Sec. 1 1 1 1 1 we present the CP- 
violating invariant products of coupling constants which 
enter into the primordial lepton asymmetry, while, in 
Sec. IIVI we sketch the calculation of the BAU. The ef- 
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TABLE I: U(l) charges of superfields 
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fccts of R-symmetry violation are discussed in Sec. [V] 
The constraints from neutrino masses and the preserva- 
tion of the initial lepton asymmetry are given in Sec. lVIl 
and our numerical results in Sec. IVIII Finally, our con- 
clusions are summarized in Sec. IVIIII 

II. THE MODEL 

The model has gauge group SU(3) C x SU(2) t xU(l)yX 
U(l) b-l and a global R-symmetry U(1),r, which is 
though explicitly broken by some terms in the superpo- 
tential (see below) . In addition to the corresponding vec- 
tor superfields and the usual MSSM chiral superfields hi, 
h 2 (Higgs SU(2) L doublets), k (SU(2) L doublet leptons), 
e\ (SU(2)x, singlet charged leptons), qi (SU(2)i doublet 
quarks), uf, d\ (SU(2)x, singlet quarks) with i = 1,2,3 
being the family index, we introduce chiral superfields 
tf (RHNs), S, (f), ({> singlets under the SM gauge group 
and T, T in the adjoint representation of SU(2)i with 
Y = 1, — 1 respectively. (As usual, we will use the same 
symbol to denote the superfield and its scalar component, 
the distinction being clear in context.) The charges un- 
der U(1)b_l and U(1)h are given in Table I. 

The superpotential is 

W = M 2 ) + \S(hih 2 ) 

+ h ei j(hili)ej + h u ij(h 2 qi)uj + h dij (hiqi)dj 

+ h„ij(h2k)Vj + hnjheTlj + h T h 2 eTh 2 

+ [M„h j /M 2 ) 2 vfv? + {M T /M 2 )4> 2 Tf 

+ (X t /M s )^hih 2 X + (\' t /M s )0hieTk + • • • , (1) 

where M is a mass parameter of the order of the GUT 
scale, Ms is the string scale ~ 5 x 10 17 GeV, e is the 2x2 
antisymmetric matrix with ei 2 = 1, (XY) indicates the 
SU(2)i invariant product e a bX a Yi, and T = T a a a /V2, 
T = T a a a j\f2 with a a (a = 1, 2, 3) being the Pauli ma- 
trices. The ellipsis represents terms of order higher than 
four and summation over indices is implied. The only 
U(1)b violating terms which we allow in the superpoten- 
tial are the two explicitly displayed terms in the last line 
of the right hand side (RHS) of Eq. QJ, which are neces- 
sary for leptogenesis. We can show that baryon number 
(B) is automatically conserved to all orders as a conse- 
quence of U(1)_r. The argument goes as in Ref. and 
is not affected by the presence of the above U(l)# break- 
ing superpotential terms. Lepton number (L) is then also 
conserved as implied by the presence of U(1)b-l- 

The inflationary trajectory is as described in Ref. 0: 
for k < A, it is parametrized by S, \S\ > S c — M, with 



all the other fields vanishing and has a constant energy 
density k 2 M 4 at tree level. Here, the dimensionless pa- 
rameters k, A and the mass M are taken real and pos- 
itive by redefining the phases of the superfields. There 
are radiative corrections [19J which lift the flatness of this 
classically flat direction leading to slow-roll inflation un- 
til \S\ reaches the instability point at|5| = M as one can 
deduce from the e and rj criteria |25j. The quadrupole 
anisotropy of the CMBR and the number of e-foldings 
N Q ~ ln[1.88 x 10 11 K 1 / 3 (A//GeV) 2 / 3 (T rh /GeV) 1 / 3 ] H| 
of our present horizon scale during inflation are given by 
the Eqs. (2)- (4) of Ref. j2(| with the two last terms in the 
RHS of Eq. (3) divided by two since the SU(2) fi doublet 
chiral superfields l c , l c of this reference are now replaced 
by the SM singlets 4>, 4> (see also Ref. p^l. 

When the value of |5| falls below M, a B — L breaking 
phase transition occurs provided that n < A. The fields 
evolve towards the realistic SUSY minimum at (S) = 0, 
((f>) = (</») = M, (hi) = (h 2 ) = 0, where ((f)), (<f>) are taken 
real and positive by a B — L rotation (there is also an un- 
realistic SUSY minimum which is given below) . Actually, 
with the addition of soft SUSY-breaking terms, the posi- 
tion of the vacuum shifts Hj to nonzero (S) ~ —^13/2/^5 
where to 3 / 2 is the gravitino mass, and a small effective 
/j, term with /1 ~ —\m 3 / 2 / K is generated from the super- 
potential coupling \S(hih 2 ). Subsequently, the inflaton 
degrees of freedom S and 6 = (5(f> + 5(f))/V2 (S(f> — cf)—M, 
8<j) = <f> — M) with mass minf = \/2kM oscillate about 
this minimum and decay to MSSM degrees of freedom re- 
heating the universe. The predominant decay channels of 
S and 6 are to fcrmionic and bosonic hi, h 2 respectively 
via tree-level couplings derived from the superpotential 
terms XS(hih 2 ) and nScfocf). Note that, if k > X, the sys- 
tem would end up in the unrealistic SUSY minimum at 
(f) = = 0, \hi\ = \h 2 \ — (k/A) 1 / 2 M, which is degenerate 
with the realistic one (up to m^, 2 ) and is separated from 
it by a potential barrier of order m^^M 2 . 

The RHNs and the SU(2)x triplets acquire masses 
Mi/cjj and Mt respectively after the spontaneous break- 
ing of U(1)_b_l by (4>), (<p). The terms which appear in 
the fourth line of the RHS of Eq. can also be written in 
the form \„ Hj 4> 2 vfvpMs, X T (f> 2 TT/M s making it clear 
that the RHN and SU(2)l triplet masses are suppressed 
by a factor MjMg relative to M. It is possible to rede- 
fine superfields to obtain effective mass terms M„ci/?i/9 
(which are diagonal in the flavor space) and MtTT with 
M„c and Mt real and positive. 

Similarly, after the U(1)b-l breaking, the explicitly 
displayed terms in the last line of the RHS of Eq. (JTJ 
which violate U(l)# give rise to effective B — L and mat- 
ter parity violating operators ^(/ii/ia)^ arid QhieTli, 
where Q and £i are suppressed by one power of M/Ms- 
If we require that the magnitude of the dimensionless 
coupling constants Xi and X[ is less than unity, we obtain 
the bound \&\, \Q < M/M s . Note that, although Q 
can be made real and positive by redefining the phase of 
li, there is no phase freedom left which can do the same 
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for Ci- This can be shown by considering the rephasing 
invariant Q 2 fj, 2 M u c (no summation over i) with /i and 
M„c already made real. The coupling constants Ci, thus, 
remain in general complex. It is, of course, possible to 
write down many other R-symmetry violating operators. 
However, they are not necessary for a nonzero lepton 
asymmetry to be created. 

After electroweak symmetry breaking, the triplet T 
also acquires HQ a VEV (T) = -hf (h 2 ) 2 /M T . This is 
due to the mass term M T TT and the coupling hfh 2 eTh 2 , 
where hj< is made real and positive by a redefinition of 
the phase of T and a compensating rephasing of T and 
li in order to retain the positivity of Mt and £j. Note 
that (hi), (h 2 ) can be taken real because of the reality of 
B[i ~ — 2Xm^, 2 /K Q. So (T) is also real. 

Performing appropriate flavor rotations of the li and 
e°p the Yukawa coupling constant matrix h ei j can be di- 
agonalized with real and positive entries h e i in the di- 
agonal. The neutrino components (z/j) of U are then in 
the weak interaction basis. (The rephasing of li which 
was used to make Q real and positive should actually 
be performed after these rotations and should be accom- 
panied by a compensating rephasing of so that h e i 
remains real and positive.) The coupling of T to a pair 
of lepton SU(2)l doublets cannot be made diagonal since 
no further flavor rotations of the k are allowed. More- 
over, if we define superfields such that Mt, Cl, A* and 
are real then it is in general not possible to make the 
h-Tij real since, once S, 0, (f> acquire nonzero VEVs, we 
have rephasing invariants MThTijh^ii 2 ^ £j* (no sum- 
mation). Finally, the Yukawa coupling constants h v ij 
also remain in general complex as one can easily deduce 
from the rephasing invariants h*p : £-* 2 M T y? M^h 2 ^ (no 

summation over indices). 

The calculation of lepton asymmetry produced in S 
and decays is quite straightforward but somewhat 
lengthy, and differs in detail from the usual case where 
leptogenesis occurs via the decay of RHN and/or SU(2)z, 
triplet superfields. Since we consider the interference of 
two one-loop diagrams, we will need to calculate the real 
parts of loop integrals which require renormalization. 

III. CP- VIOLATING INVARIANTS 

To produce a nonvanishing net B — L asymmetry from 
inflaton decay, the theory must contain one or more phys- 
ical CP-violating quantities. These are products of cou- 
pling constants, corresponding to operators noninvariant 
under CP conjugation, which are nonreal and are not af- 
fected by the redefinition of fields by complex phases (or 
other global symmetries). In this case, we consider the 
terms of W in Eq. . Since leptogenesis takes place at 
reheating, we work in the vacuum where (</>) = (4>) = M. 
For simplicity, we ignore the couplings of the inflaton to 
RHNs and SU(2)l triplets from the terms in the fourth 
line of the RHS of Eq. JTJ together with its couplings re- 
sulting from the two subsequent terms. The inclusion of 



these couplings would only complicate the analysis with- 
out altering the character of our mechanism. So these 
four superpotential terms are replaced by the effective 
mass terms M^ufuf, M^TT and couplings C,i{h\h 2 )v^ 
and QhicTli. The condition of CP invariance is simply 
that it should be possible to make all products of cou- 
pling constants real by field redefinitions. 

We consider here only rephasing invariants which, in 
contrast to the standard invariant used in Ref. Q , would 
exist even if there was a single fermion family. We find 
the following independent invariants (no summation over 
repeated indices): 

ha = M T huijCj{Mu9hfCi)*, 
hijk — MThvijhvkjiMvcfiTikhf)* , (2) 

whose imaginary parts violate CP invariance. Any com- 
bination of the two (modulo real invariants) results in a 
third invariant such as I^ijk = C'^fcj (Cj^Tifc)* ■ In each 
case, the expression involves both the SU(2)i triplet and 
RHN couplings, thus the generation of a B — L asym- 
metry is independent of the sources of CP violation con- 
sidered in previous scenarios and we require novel decay 
diagrams. The invariants Iuj, I 2 ijk, l3ijk are the mini- 
mal ones in the sense that they have the smallest possible 
number of trilinear superpotential couplings. 

It is important to note that the invariant Iuj can be 

split in two parts C*Mxh%, and QM*ah V ij correspond- 

j 

ing to effective operators which carry (opposite) nonzero 
B — L charges and involve only light fields since the heavy 
ones can be contracted. These light fields include bosonic 
or fermionic hi, h 2 and their conjugates. So, Iuj is, 
in principle, suitable for leptogenesis which requires the 
interference of two B — L violating diagrams. (Recall 
that the inflaton field couples at tree level to the elec- 
troweak Higgs superfields hi, h 2 and can decay only to 
light particles.) This is not the case with the rephasing 
invariant I 2 ijk since there are no hi's among the light 
fields of the corresponding (effective) operator. On the 
contrary, the invariant I^ijk has all the above good prop- 
erties of Iuj , but one of the interfering amplitudes turns 
out to be zero in this case. The reason is that it in- 
volves a real (with vanishing absorptive part) on-mass- 
shell off-diagonal self-energy between / and hi which, in 
the on-shell (OS) renormalization scheme, vanishes (sec 
Sec. lIVII . It is clear that any invariant which can be useful 
for leptogenesis must involve Iuj and, thus, the effective 
coupling constants Q, So, the explicit violation of 
matter parity is essential for our scheme. This is another 
novel feature of this leptogenesis scenario. 

IV. BARYON ASYMMETRY 

The CP-violating rephasing invariant Iuj corresponds 
to the product of coupling constants in the interference 
of the diagrams in Fig. QJi (Fig. with the diagrams 
in Figs. ^) and^ (Figs. [2b and|2t) for the L violating 
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FIG. 1: The nine one-loop diagrams for the L violating decay 
S — > U /i2- The solid (dashed) lines represent the fermionic 
(bosonic) component of the indicated superfield. The arrows 
depict the chirality of the superfields and the crosses are mass 
insertions in fermion lines. 



decay of S (9). This interference contributes to the L 
asymmetry due to a partial rate difference in the decays 
S -> l t h 2 and S* -> I* h* (9* -> k~h 2 and -> khv), 
where bar and tilde represent the antifermion and the 
SUSY partner respectively. Both the decaying infiaton 
(S or 9), which we take at rest, and the decay products 
must be on mass shell. For simplicity, we consider that 
all the propagating and external MSSM particles in the 
diagrams are massless. Also, we perform the calculation 
in the limit of exact SUSY. 

In each case, the resulting contribution to the L asym- 
metry is proportional to both Im Iuj and the imaginary 
part of the interference of the relevant 'stripped' dia- 
grams with the dimensionless coupling constants and the 
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FIG. 2: The three one- loop diagrams for the L violating 
decay 9* — > Uhz- The notation is as in Fig. 



My, M u c mass insertions factored out (we keep, though, 
the minf factor appearing in the scalar coupling 9*h\h 2 ). 
The 'stripped' diagrams in Figs. HJi, ^ and^; are de- 
noted by Ff m , Fg jn and F% jn respectively with i and j 
being the family indices and n = 1,2,3 the serial num- 
ber of the diagram. Similarly, the 'stripped' diagrams 



in Figs. Hi, Hd and|2i are F^, F h 
total net L asymmetries e\g and e 
9 decay respectively are given by 



bij 



and F e clJ 



Thus, the 
generated per S and 



e |s - -2^-Im/iijIm [F ain F U j n +F ain F cijn ], 

e\ e = -2^lmIuM[KiKi* + KiFtijX (3) 

where T = |A| 2 mi n f/87r is the rate of the tree-level decays 
S — > h\h 2 and 9 — > h\h 2 , and summation over the in- 
dices i, j, n and integration over the phase space of the 
final particles is implied. 

We see that the one-loop diagrams for the L violat- 
ing decay of the infiaton which contain SU(2)j, triplet 
exchange are exclusively of the vertex type (see Figs. ^ 
and[5Ji). On the contrary, the diagrams with a RHN 
exchange are both of the vertex (Figs. QJj and [2b) and 
self-energy (Figs.^J; and|2t) t 2§] type. Each of the three 
vertex diagrams in Fig. ^ (Fig. QJ)) possesses a logarith- 
mic ultraviolet (UV) divergence. However, one can easily 
show that their sum equals TOi n f times the vertex diagram 
in Fig. [2]} (Fig. Ht>)> which is UV finite. Similarly, one 
can show that, although each of the three self-energy di- 
agrams in Fig. 2J: has a quadratic divergence, their sum 
equals the self-energy diagram in Fig. [2t multiplied by 
minf. However, the latter is not UV finite. It rather pos- 
sesses a logarithmic divergence and, thus, needs renor- 
malization (see below). 

The above relations between the diagrams for the L 
violating decays of S and 9 imply that e\$ = e\e = £• 
We can, thus, concentrate on the calculation of e\$ which 
is simpler. As already mentioned, the vertex diagrams 
in Figs. and \Bp are finite (both their real and imag- 
inary parts) and, thus, well-defined and independent of 
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the renormalization scheme used. However, the diagram 
in Fig. involving a divergent self-energy loop requires 
us to apply a renormalization condition. As shown in 
Ref. the appropriate renormalization scheme, in this 
case, is the on-shell (OS) scheme. In a general theory 
with scalars Si, the OS conditions on the renormalizcd 
self-energies ILjj(p 2 ) are as follows: 



ReU lJ { f if)=ReU lJ (fij) = 
for the off-diagonal self-energies (i j) and 
1 



lim 



Ml 



rRelLuip- 1 ) = 



(4) 



(5) 



for the diagonal ones (see e.g. Ref. 29]). Here we take 
a basis where the renormalized mass matrix is diago- 
nal with eigenvalues /i^. The imaginary part of the self- 
energies is finite and, thus, not renormalized. In Fig. [2t, 
we have an off-diagonal self-energy diagram between the 
scalars 9 and z>|. Given that 9 is on mass shell, the real 
part of this diagram vanishes in the OS scheme. The 
imaginary part, however, gives a finite contribution. 

The lepton asymmetry per inflaton decay is calculated 
by making use of the software packages of Ref. [3(j and 
is found to be given by 



evv + evs, 



where 



Im/ 



evv 



Im [f(yr)f(yj 



(6) 



(7) 



"inf 



is the contribution from the interference of the two vertex 
diagrams in Figs. and [2t> , and 



Im / 



CVS 



64tt 3 m? n{ - Ml 



Re/(?/7 



(8) 



is the contribution from the interference of the ver- 
tex and self-energy diagrams in Figs. and with 
f(y) = n 2 /6-U 2 (l + y + ie) (Li 2 is the dilogarithm 31]), 
Ut = mLt/Mq,, yj — m? JM^ and summation over the 

family indices i, j implied. It should be emphasized that 
Eq. JSJ holds [TjJ provided that the decay width of is 
<C \mf n[ — Mlc\/nii n f, which is well satisfied in our model 

3 

if M„= is not unnaturally close to mi n f. 

The equilibrium conditions including nonpcrturbative 
electroweak reactions above the critical temperature of 
the electroweak phase transition yield a relation between 
the baryon number density ub and the B — L number 
density ub-l, which allows us to find ub in terms of the 
tib—l produced from inflaton decay, assuming that B — L 
is conserved at temperatures well below My and M„<= (see 
Sec.|Vl|>- In the MSSM with soft SUSY-breaking terms, 
we have tib/s = (28/79)ns_i/s 32]. If we imagine the 
inflaton to decay instantaneously out of equilibrium cre- 
ating initial lepton number density 71l, init then 



s 



28 rt£,init 
79 s 



28 Uin{_ 

79 £ s 



0) 

79 minf 



using the standard relation n- ul {/s = (ns + ng)/s — 
3T r h/4mj n f for the inflaton number density. Here s is 
the entropy density. The reheat temperature is given by 



T, 



rh 



/ 45 



(Tmp)' 



(10) 



where g* counts the relativistic degrees of freedom taking 
account of the spin and statistics and is equal to 228.75 
for the MSSM spectrum. 



V. EFFECTS OF R-SYMMETRY VIOLATION 

One would expect the explicit violation of U(l)# in the 
superpotential to have important consequences for phe- 
nomenology and cosmology. The Z2 subgroup of U(l)/j, 
which is left unbroken by the soft SUSY-breaking terms, 
is called matter parity since all the matter (quark and 
lepton) superfields change sign under it. Combined with 
the Zi fermion parity (under which all fermions change 
sign) yields R-parity, which, if unbroken, guarantees the 
stability of the LSP. In our model, however, matter parity 
is violated along with the U(l)# by the two explicitly dis- 
played terms in the last line of the RHS of Eq. Q , which 
are needed for generating an acceptable BAU. Thus, R- 
parity is broken and, consequently, the LSP generically 
becomes unstable and decays rapidly, rendering it unsuit- 
able for the role of dark matter and leading to distinctive 
collider signatures. 

In our case, if the LSP contains a Higgsino component, 
it could decay into a pair of electroweak Higgs bosons 
and a lepton. The dominant diagrams are constructed 
from the U(1)a and R-parity violating Yukawa vertices 
Q{hih2)v'j or ^hieTli with the fermionic Vj or T con- 
nected to the fermionic or T of the Yukawa couplings 
h l iij{h,2li)v'j or h'fhitThi respectively via a mass inser- 
tion. For the numerical values of the parameters that we 
consider (see Secs. lVll and lVll|l . we find that the LSP life- 
time can be as low as about 10 _1 sec. However, it is easy 
to block kinematically the LSP decay by taking its mass 
to be smaller than twice the mass of the lightest Higgs 
boson, which is very reasonable. Thus, it is perfectly 
possible to rescue the LSP as dark matter candidate. 

Our model could also predict the existence of other 
R-parity violating processes at low energies besides the 
LSP decay. The R-parity breaking superpotential cou- 
plings involve at least one superheavy field (a RHN or 
SU(2)l triplet). On integrating out these heavy fields, 
one generally obtains effective R-parity violating opera- 
tors involving only MSSM fields. These operators, if they 
have dimension five or higher, do not lead to detectable 
processes since they are suppressed by some powers of 
Mt or M„<:. However, dimension four operators such as 

the effective scalar vertices /ii/i^/i^i* or hikljlk' which 
originate from the superpotential couplings Q(hih2)vj, 
hvij(h2h)vj or ^hieTli, fiTjkljcTlk respectively, can lead 
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to low-energy R-parity violating processes which may be 
detectable in the future colliders (see Sec. IVll| . 

It is well known [33| that, in any leptogenesis scenario 
with RHNs or SU(2)l triplets, it is important to ensure 
that the primordial lepton asymmetry is not erased by 
lcpton number violating 2^2 scattering processes such 

as lilj — > h 2 h 2 or lih 2 — > h 2 l* at all temperatures be- 
tween T r h and about 100 GeV. In our model, due to 
the presence of the R-parity violating superpotential cou- 
plings, there exist some extra processes of this type such 



as h\L — > h 2 h 2 or h 1 h 2 



h*l* 



which are derived from 



diagrams similar to the ones mentioned above for the LSP 
decay. In addition to all these processes which correspond 
to effective operators of dimension five (or higher), we 
also have dimension four R-parity (and lepton number) 
violating processes such as h\h 2 — > h 2 h or h\U — > Ijlk, 
which are derived from the effective four-scalar vertices 
in the previous paragraph. 

The initial lepton asymmetry is protected |32j| by 
SUSY at temperatures between T r h and about 10 7 GeV. 
For T < 10 7 GeV, one can show that all the lepton num- 
ber violating 2^2 scattering processes which result 
from effective operators of dimension five or higher are 
well out of equilibrium for the values of the parameters 
used here. The dimension four lepton number violat- 
ing processes, however, are generally in equilibrium and, 
thus, special care is needed in order to retain the initial 
lepton asymmetry in our model. We will return to this 
issue in the next section. 

As already explained, the classical flatness of the infla- 
tionary path in the limit of global SUSY is ensured, in our 
model, by a continuous R-symmetry enforcing a linear 
dependence of the superpotential on S. This is retained 
[22| even after SUGRA corrections, given a reasonable as- 
sumption about the Kahler potential. The solution to 
the /i problem is also reliant on the R-symmetry. These 
aspects of the model are not affected by the explicit R- 
symmetry breaking we consider. 

In the model, some R-symmctry violating couplings 
are present in the superpotential and some not. Thanks 
to the nonrenormalization property of SUSY, this situ- 
ation is stable under radiative corrections, but one may 
consider it unnatural since there is no symmetry to forbid 
the terms we set to zero. 



VI. NEUTRINO MASSES AND PRESERVATION 
OF LEPTON ASYMMETRY 

The two distinct sources of neutrino masses in the 
model yield the following mass matrix of light neutrinos: 



= -(h 2 



hTi i i h 1 h 



hrpMx 



(11) 



The terms in this expression are generally complex with 
physical relative phases between the two contributions. 



Comparing this to the combination of parameters ap- 
pearing in the lepton asymmetry per inflaton decay, we 
see that there is no obvious correlation between the two. 
Hence, we have here more freedom to choose values for 
the parameters than in the case of leptogenesis with a 
single source of neutrino mass. However, this freedom is 
reduced by requiring that the lepton asymmetry is pre- 
served until the electroweak phase transition. 

As discussed in Sec. there exist unsuppressed di- 
mension four R-parity and lepton number violating 2 — ► 2 
scattering processes which can erase the primordial lep- 
ton asymmetry. They originate from the effective four- 
scalar operators h\h 2 h 2 l* or hxliVjVf. with coupling con- 
stants J2j Cjhtij or Cih-Tjk respectively. The former oper- 
ator violates, in particular, the i-lepton number (Li), i.e. 
the part of the total lepton number which corresponds 
to the ith fermion family. Therefore, in order to retain 
at least one part of the lepton asymmetry, we need to 
impose the condition that 



h* ■ 

3 11 vi j 



= 0, 



(12) 



for some value of i [34| . To generate a Li asymmetry, we 
need to take a nonzero £[, as it is obvious from Eqs. (J2J 
and The operators hxkf-1% with j,k ^ i or j = 

k = i will then violate Li, thereby leading to the erasure 
of the primordial Li asymmetry. To avoid this, we must 
take IiTjk to vanish, unless exactly one of the indices is 
equal to i, i.e. we allow only hxik = hxki ^ for k =/= i. 
Moreover, if fiTik tnen we require £j = for j : ^ i to 
avoid generating the Li violating operator h\ljl*l^. 

We must further ensure that there are no processes in 
equilibrium which violate Li, while conserving L. Such 
processes result from the four-scalar operators h 2 lkh 2 l* 
(k =/= i) which originate from the superpotential coupling 
h v kj{h 2 lk)i / j ■ Thus, we must take 



^ Kkjh* 



3 



0. 



(13) 



for the value of i which satisfies Eq. I|12f) and all values 
of k which are different from this i. Also, the operators 
hilih^lj with j ^ i, which originate from the superpo- 
tential coupling (jhieTlj, violate L i: but not L. So, we 

must take Cj — for j ^ i. Finally, the operators Ijlkl^n 
derived from the superpotential term hxjklj^Tlk could 
violate Li. However, as explained above, hxjk should 
vanish unless exactly one index is equal to i. Thus, the 
remaining operators automatically conserve Li. In sum- 
mary, for an asymmetry to be generated in Li and not 
wiped out by scattering processes, we require that ^ 
and = for all j ^ i. Also, hrjk = unless exactly 
one of j or k is equal to i and Eqs. (|12|) and 113|) should 
hold for all k ^ i. 

We take the neutrino mass ordering mi < m 2 < 
and adopt the normal hierarchical scheme of neutrino 
masses, where the solar and atmospheric neutrino mass 2 
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differences are identified with Srn^i and 6m 2 .-, respec- extra conditions: 
tively. Analysis [35| of the CHOOZ experiment |36j shows 
that the mixing angle #13 can be taken equal to zero. 
Moreover, the fact that &m\ x <C 8m\^ implies that, when 
considering atmospheric neutrino oscillations, it is a good 
approximation to set bm\ x = 0. For simplicity, we fur- 
ther take a = (3, where a, (3 are the Majorana phases 
in the leptonic mixing matrix associated with the two 
lighter neutrino mass eigenstates P±, vi respectively. 

Under these circumstances, the weak interaction eigen- 
state v\ decouples from the other two ^2, ^3 in the neu- 
trino mass matrix. A simple choice of parameters which 
leads to this decoupling is h v ij = for j = 2, 3, h u n = 
for i = 2, 3 and hxn = hxij — h-Tji = for j — 2, 3. 
Eq. 1)13(1 is then automatically satisfied for k = 1 and 
i = 2,3, which means that processes corresponding to 
the operators h^lxh^l* (i — 2,3) which violate L\, while 
conserving L, are avoided. If we further take Ci = 0; 
L\ violating and L conserving processes from the oper- 
ators hihhllj (j — 2,3) are also absent. Moreover, our 
choice of fiTjk ensures that no such processes from oper- 
ators of the type Ijlklm^n appear. Eq. (|12fl with i = 1 
(i = 2,3) receives contribution only from the term(s) 
with j = 1 (j — 2,3). Also, L and L\ violating pro- 
cesses from the operators hilil*-l%. do not exist for our 
choice of parameters. Note that there are no L\ vio- 
lating processes in equilibrium from the superpotential 
terms ^hieTl^, hxjklj^Tlk- The only possible L\ vio- 
lating processes come from the operator hihsh^l*, which 
is though L2 and L3 conserving. The fact that C,[ = 
implies that no primordial L\ asymmetry is generated as 
one can see from Eqs. @ and Also, the L2 and 

L3 asymmetries cannot turn into L\ asymmetry by 2 — > 2 
scattering processes since, as explained, there are no such 
processes in equilibrium which simultaneously violate L\ 
and L2 or L3. All these facts allow us to ignore the first 
family and concentrate on the two heaviest families. 

From the recent global analysis [37j of neutrino oscil- 
lation data, we take the best-fit value of the atmospheric 
neutrino mass 2 difference which, in the case of the normal 
hierarchical scheme, yields the value 7713 ~ 5.1 x 10 -2 eV 
for the heaviest neutrino mass. Also, the mixing angle 
023 is identified with its best-fit value which is about 7r/4 
and 777,2 is taken equal to zero consistently with the fact 
that we are considering only atmospheric neutrino oscilla- 
tions. So, the neutrino mass matrix elements in Eq. 1(11(1 
must satisfy the following restrictions: 



K32 + RK33 - °> 



^22^32 



^23^33 — 0' 



(15) 



where R — C,z/C,2- Applying the results of the previous 
discussion, we require £3 to be nonzero and C2 = (so 
no L2 asymmetry is generated); we have also hx22 = 
hx33 = 0, so the only nonzero coupling constant of the 
type fiTjk is hx23 = ^T32- Given the value of the right 
handed neutrino masses, Mt and hf, the five conditions 
in Eqs. ((14(1 and 1(15(1 can be solved, using the mass for- 
mula in Eq. 1(11(1 . to find a relation between the complex 
Yukawa coupling constants hij (i, j = 2, 3) and the com- 
plex parameters R and hx23- 

Note that the presence of L and B — L violating 2-^2 
scattering processes leads to a moderate modification of 
the numerical factors in Eq. ©. However, this modifica- 
tion depends on details and we will, thus, ignore it. 



VII. NUMERICAL RESULTS 

We saturate the gravitino bound on the reheat temper- 
ature by taking T r h = 10 10 GeV which is acceptable 01 
provided that the branching ratio of the gravitino decay 
to photons and photinos is less than unity. We also fix the 
parameter k to the value 10 -4 . The cosmic microwave 
background explorer (COBE) value of the quadrupole 
anisotropy of the CMBR ((5T/T) Q ~ 6.6 x 1CT 6 ) [H 
is then reproduced for A ~ 1.44 x 10~ 4 (> K as it should) 
and M ~ 3.53 x 10 15 GeV. Thus, 777 inf ~ 4.99 x 10 11 GeV. 
The spectral index of density perturbations comes out 
practically equal to unity. 

We now wish to demonstrate that our leptogenesis 
mechanism can easily reproduce the best-fit value of the 
BAU, tib/s ~ 8.66 x 10~ n , derived from the recent 
WMAP data Q with a natural choice of parameters. For 
simplicity, we approximate all the right handed neutrino 
masses to the same value M„c. (Note that we do not 
make use of resonance effects for degenerate heavy neu- 
trinos, so this approximation will not strongly influence 
our results.) Eq. I(15|) implies that h V 32 = —R*h v ^ and 
h V 23 — Rh U 22- (Note that /i„33 7^ for a nonvanishing 
L3 asymmetry to be produced.) Substituting these in 
Eq. I(14|). we then obtain 



1 

WV22 = «V 33 = "V23 = 2 m 3- 



±i 



(14) 



7773 M v c 
2(/7 2 } 2 J (l + i? 2 )i 



(16) 



(Note that, in the case of three neutrino flavors, our 
choice of parameters supplemented, for consistency, with 
a vanishing mi yields the extra restriction h v \\ = 0.) 
These equalities need not be exact, however, we will take 
them as so for simplicity. We choose to maintain the L 3 
asymmetry, which requires that Eq. ((12(1 holds for i = 3 
and Eq. ((13(1 for k = 2, i = 3. So, we obtain the following 



and 



± 



J i/33 



±i 



R-R 



{1 + R 2 )i(l + R* 2 )i 



\2(h 2 ) 2 J (l + i?*2)i 

2/7 T23 (/72) 2 

hf 7773 Mt 



(17) 



(18) 



9 



Note that the signs of h U 22 and can be chosen inde- 
pendently; however, the sign on the left hand side (LHS) 
of Eq. (|18|l is the product of the two. 

The LHS of Eq. I|18|) is imaginary. Therefore, the real 
part of the RHS of this equation must vanish, which gives 



100 



Re hm = - 



hrpm^MT 
2(h 2 ) 2 ' 



(19) 



This implies that 
R — R* 



(l+i? 2 )2(l+i?*2) 3 

which yields 

£ 2 |i?| 4 + 2[(1 + C 2 ) cos2<^ - l]|i?| 2 + i 2 = 0, (21) 

where tp is the phase of the complex parameter R. This 
equation admits real solutions for |i?| 2 provided that 
-1 < cos2<yS < (1 - £ 2 )/(l + £ 2 ). They turn out to 
be nonnegative and are given by 

ii?i 2 = r 2 {[i-(i+acos2^]±(i+a^ 

[(1 + C 2 ) cos 2 2<p - 2 cos 2ip + (1 - £ 2 )] 5}. (22) 

From Eqs. © and ©-©, we see that, in our case, the 
baryon asymmetry is proportional to C3^TlC2^33l m R\ 
for given M v o,M T > m inf /2 ~ 2.5 x 10 11 GeV. Using 
Eqs. (|17|) and l)2()[l. we can further show that this expres- 
sion is proportional to ^hf^S, 1 ^ 2 R x ^ 2 \i^ — cos2ip) 1 / i . 
As discussed in Sec. [H] the parameters ('3, IC2I and 
|Cs| = l-HUCal should not exceed M/M s ~ 7 x 10~ 3 . To 
maximize ng/s, we saturate this limit on £3 and on |£ 2 | 
or |^3 1 for < 1 or \R\ > 1 respectively. 

In the former case, the baryon asymmetry, which is 
oc hf\£\ 1 / 2 \R\ 1 / 2 (l—cos2(p) 1 / 4 ', can be further maximized 
by choosing the positive sign in the RHS of Eq. H22J1 and 
taking cos2<p = -1. We then get \R\ = (1 + S,- 2 ) 1 / 2 + 
which is always > 1. Our requirement that \R\ 
does not exceed unity is fulfilled only asymptotically, i.e. 
for |£| — > 00, where \R\ — > 1. The baryon asymmetry 
is then cx hf |£| 1 / 2 and is maximized by maximizing h T 
and |Im/iT23|- It would thus be desirable to put both 
these quantities equal to unity. The parameter hf can be 
readily fixed to unity and this improves the naturalness 
of our scheme. However, if |Imft,T23| is too large, some 
moderate cancellations between different contributions to 
the neutrino mass matrix are required. So, we take values 
of |Im/iT23| which are smaller than unity, but greater 
than, say, 10~ 2 to be consistent with the requirement 
that |£| ^> 1 for Mt not much bigger than mj n i/2. 

In the case where \R\ > 1, we should saturate the limit 
on I £3 1 rather than the limit on IC2I, which yields that the 
baryon asymmetry is oc /iy|^| 1 ' 2 |i2| _1 ' 2 (l — cos 2<^) 1 / 4 . 
This is maximized by choosing the negative sign in the 
RHS of Eq. IL'l't and taking again cos2(,c = —1, which 
gives \R\ = (1 -f^" 2 ) 1 / 2 — |£| _1 < 1. So, marginal consis- 
tency is achieved again for |£| ^> 1, where \R\ — > 1. The 




4 5 6 7 8 9 

M v c (10 11 GeV) 



10 11 



FIG. 3: The solutions of n B /s = 8.66 xltT 11 in the M^-M T 
plane for Im/iT23 = 0.02 (solid lines) or 0.04 (dashed lines). 
The values of all the other parameters are given in the text. 



baryon asymmetry is again oc /i-r|^| 1/ ' 2 and is maximized 
by taking hj. = |Im/iT23| = 1- As in the previous case, 
we take hf = 1 and |Im/iT23| ~ 10~ 2 . The result is ob- 
viously the same in both cases to a good approximation. 
The phase of £2 is appropriately adjusted in each case so 
that we obtain the maximal baryon asymmetry. 

For each M v c, Mt > m- m f/2, we calculate the BAU us- 
ing Eqs. I©-© with Im/i T 23 = 0.02 or 0.04 and values 
for the other parameters as explained above. We then 
compare the result with the best-fit value of the BAU 
from WMAP. The resulting solutions in the M„c — My 
plane are shown in Fig. |3| with the solid or dashed lines 
corresponding to Im/iT23 = 0.02 or 0.04 respectively. 
Note that there are two branches in each case, one with 
A/„c > mi n f and one with M v c < m; n f. They correspond 
to the two possible signs of the denominator in the RHS 
of Eq. © for eys- Actually, this quantity, as it involves a 
s-channel exchange of a v c boson, can be easily enhanced 
by letting mj n f approach the j/ c -pole. This fact assists us 
to achieve the WMAP value of ub / s with natural values 
of the parameters. It is important though to stress that 
there is no need to get too near the ^ c -pole for reasonable 
values of Mt- In fact, the solution is very far from being 
unnaturally close to this pole. 

As mentioned in Sec. the explicit R-parity viola- 
tion in our model, which is required for leptogenesis, has 
some low-energy signatures which come from dimension 
four effective scalar vertices and may be observable in 
the future colliders. Such signatures may typically be 
the three-body slepton decay processes 



h 



hxh^hX and I2 — * hil^l 



3 f 3' 



(23) 
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which can easily be kinematically allowed. (Note that, 
for our choice of parameters, similar l\ decay processes 
do not appear, although they may be present in the gen- 
eral case.) The effective coupling constants of the pro- 
cesses in Eq. iffify are + CzKs23 an d Ci^to re ~ 
spectively. Using the relation h U 23 = Rh V 22 and sub- 
stituting h V 22 from Eq. we can easily show that 
the magnitude of the former coupling constant becomes 
~ |C2||2Im/i T23 | 1 / 2 (M ly c/M T ) 1 / 2 in the large |£| limit. 
For Im/ir23 = 0.02 and taking M v c ~ Mr, which holds 
near the right corner of the solid line in Fig. |3J we find 
that the magnitude of the effective coupling constants of 
the processes in Eq. is about 1.4x 10 -3 and 1.4x 10~ 4 
respectively. The corresponding decay rates are then of 
order 10 -8 GeV and 10~ 10 GeV respectively for mass of 
the decaying slepton ~ 1 TeV and assuming that there 
is an appreciable gap between this mass and the sum of 
the masses of the decay products. 

It is finally interesting to point out that, as one can 
readily show, our scheme fails to provide any useful pre- 
dictions for the so far undetermined or not so accurately 
determined parameters in the neutrino mixing matrix, 
i.e. the three CP-violating phases and the mixing an- 
gle 0i3 . The reason is that the number of parameters is 
such that successful leptogenesis is possible whatever the 
values of the complex phases and #13. 

VIII. CONCLUSIONS 

We proposed a scenario of nonthermal leptogenesis fol- 
lowing supersymmetric hybrid inflation, in the case where 
the decay of the inflaton to both heavy right handed 
neutrino and SU(2)l triplet superfields is kinematically 
blocked. The primordial lepton asymmetry is generated 
through the direct decay of the inflaton into light parti- 
cles. We implemented our scenario in the context of a 
simple SUSY GUT model which incorporates the stan- 
dard version of SUSY hybrid inflation. The \x problem is 
solved via a U(l) R-symmetry which forbids the existence 
of an explicit /x term, while allows a trilinear superpo- 
tential coupling of the gauge singlet inflaton to the elec- 
troweak Higgs superfields. After the spontaneous break- 



ing of the GUT gauge symmetry, this singlet inflaton ac- 
quires a suppressed VEV due to the soft SUSY-breaking 
terms. Its trilinear coupling to the Higgs superfields then 
yields a /1 term of the right magnitude. 

The main decay mode of the inflaton is to a pair of 
electroweak Higgs superfields via the same trilinear cou- 
pling. The initial lepton asymmetry is created in the sub- 
dominant decay of the inflaton to a lepton and an elec- 
troweak Higgs superfield via the interference of one-loop 
diagrams with right handed neutrino and SU(2)x triplet 
exchange respectively. The existence of these diagrams 
requires the presence of some specific superpotential cou- 
plings which explicitly violate the U(l) R-symmetry and 
R-parity. However, the broken R-parity need not have 
currently observable low-energy signatures, although it 
may have signatures detectable in future colliders. Also, 
the LSP can be made stable and, thus, be a possible 
candidate for cold dark matter. 

In our analysis, we took into account the constraints 
from neutrino masses and mixing. There exist, in our 
model, two separate contributions to the neutrino mass 
matrix, which originate from the usual seesaw mecha- 
nism and from the SU(2)l triplet superfields. The con- 
straints arising from neutrino masses and mixing alone 
are not very stringent. However, the requirement that 
the primordial lepton asymmetry not be erased by lep- 
ton number violating processes before the electroweak 
phase transition is a much more stringent constraint on 
the parameters of the theory. Taking into account these 
constraints we found that the best-fit value of the BAU 
from the recent WMAP data can be easily achieved with 
natural values of parameters. 
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